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INTRODU CTION : THE  EXPERIMENTAL  EVIDENCE 

7 

The  general  problem  of  the  interface  between  two  phases  which  are 
charged  and  or  conducting  is  of  relevance  to  a  number  of  systems  which 
occur  in  nature:  colloids,  micelles,  membranes,  solid- solution 
interefaces  in  general  and  metal  solution  interfaces  in  particular, 
form  a  bewildering  array  of  systems  of  enormous  complexity.  The 
investigation  of  the  structure  of  these  systems  poses  considerable 
difficulties,  both  experimentally  as  well  as  theoretically.  The 

O 

experimental  problem  is  that  the  interface  has  '  particles  relative 
to  the  bulk,  solid  or  liquid  phases.  For  this  reason  one  needs  a 
surface  specific  method,  which  is  able  to  discriminate  between  the 

signal  from  the  surface  from  the  rest.  Electrons  do  not  penetrate 
into  solids  and  for  that  reason  have  been  used  extensively  for  the 

determination  of  the  surface  structure  of  solids.  However  they  must  be 

used  in  vacuum  and  that  precludes  their  use  in  the  study  of  the  liquid 

solid  interface.  The  study  of  electrode  surfaces  removed  from  the 

liquid  cell  under  various  conditions  has  provided  an  enormous  wealth 

■  T " 

of  useful  data  which  we  will  not  try  to  review  here  .  The  only  way  to 

understand  the  relation  between  the  ex-situ  and  is  to  measure  both. - 

Neutrons  are  another  possible  choice:  unfortunately  there  is  no 
mechanism  by  which  the  neutrons  can  be  made  surface  selective,  and  for 
that  reason  only  systems  with  large  specific  surface,  which  in  general 
are  disordered,  can  be  studied. 

X-rays  are  a  viable  alternative  for  the  study  the  solid-liquid 
interfaces.  The  X-ray  probes  are  reviewed  in  the  excellent  chapter  by 
H.D.  Abruna.  There  are  three  techniques  that  have  been  used  to 
determine  structural  features  in  well  characterized  metal  electrolyte 
interfac  s: 

1)  EXAFS  or  extended  X-ray  absorption  fine  structure^,  which  permits 

the  determination  of  the  near  neighbor  structure  of  a  given  target 

atom,  and  also  yields  information  about  the  electronic  state  of  that 
atom  when  it  is  adsorbed  at  the  surface^. 

2)  GIXS,  or  grazing  incidence  X-ray  diffraction*.  This  method  permits 

the  determination  of  the  in-plane  structure  of  an  adsorbed  monolayer 

of  a  target  atom  on  the  surface.  It  is  a  very  exacting  technique,  and 

requires  that  the  structure  of  the  adsorbed  monolayer  should  be 


different  of  that  of  the  substrate5. 

3)  Standing  wave  methods**.  In  this  case  an  X-Ray  standing  wave  is  set 
up  at  the  interface  of  the  solid  and  the  fluid  .  There  are  several 
modes  in  which  these  standing  waves  can  be  formed;  they  allow  the 
determination  of  the  distance  from  the  surface  of  the  solid  into  the 
fluid  phased 

The  spectroscopic  methods  using  ultraviolet,  visible  or  Raman 

8 

spectroscopy  are  very  useful  in  situ  probes  .  Although  they  are  not 
directly  related  to  the  geometrical  structural  parameters,  these  can 
be  extracted  from  theoretical  considerations  with  a  fair  degree  of 
reliability.  The  optical  spectroscopic  methods  do  not  require  special 
installations  such  as  the  synchrotron,  and  are  most  useful  for  complex 

molecular  species;  The  techniques  are  the  surface  enhanced 

9  10 

Raman  ,  surface  infrared  spectroscopy  ,  second  harmonic 

generation11,  which  permits  to  discriminate  between  differer.: 

geometries  of  the  adsorbates  on  single  crystal  surfaces. 

A  technique  that  provides  with  direct  structural  information  of 

the  electrode  surface  is  the  Scanning  Tunnel  Microscope  (STM).  It  has 

been  recently  shown  by  various  groups  that  the  STM  is  capable  of 

resolving  structural  details  of  metal  surfaces  in  contact  with 

12 

electrolytic  solutions  .  However  when  the  electrochemical  potential  is 
scanned,  then  the  tunnel  voltage  of  the  STM  also  changes.  This  does 

not  affect  the  study  of  surface  geometry,  since  the  images  are 
relatively  independent  of  the  tunnel  voltage15.  The  resolution  of  the 
STM  pictures  of  the  metal  electrolyte  interface  is  of  the  order  of 
1-2  Angstroms. 

An  in  situ  electrochemical  technique  that  has  been  established 
recently  is  the  quarz  microbalance:  this  instrument  can  measure  small 

changes  in  the  mass  of  a  metallic  electrode  that  is  attached  to  a 
quarz  oscillator1^.  In  this  way  the  electrosorption  valency15  can  be 
calculated  directly  from  the  amount  of  charge  from  the,  voltamogram, 
and  the  mass  obtained  from  the  microbalance.  Proper  interpretation  of 
the  results  of  this  instrument  requires  electrode  surfaces  that  have 

large  molecularly  smooth  regions. 

A  further  method  that  has  yielded  very  interesting  information 
about  the  structure  and  interactions  in  the  diffuse  part  of  the  double 
layer  is  the  direct  measurement  of  forces  between  colloidal 

particles1**.  Here  the  forces  between  two  mica  plates  are  measured 
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directly  in  the  presence  of  different  solutions.  Quite  interestingly 

the  forces  give  raise  to  oscillations  of  a  period  similar  to  the 

dimensions  of  the  molecules  enclosed  between  the  mica  plates. 

Amongst  the  optical  techniques  there  are  also  the  more 

traditional  methods  such  as  the  ellipsometry,  electroreflectance  and 

particularly,  surface  plasmons^,  where  experimental  and  theoretical 
18 

advances  have  made  it  possible  to  offer  a  picture  of  the  surface 

electronic  states  of  the  metal  in  some  selected  cases,  such  as  the 

silver  (111)  phase.  We  should  mention  here  the  measurement  of  image 

19 

potential  induced  surface  states  by  electroreflectance  spectroscopy  . 
In  this  case,  besides  the  normal  surface  states  which  arise  from  me 

termination  of  the  crystal  lattice,  there  are  discrete  states  due  to 

the  existence  of  an  image  potential  for  charges  near  the  conducting 

interface. 

And  last,  but  certainly  not  least,  there  is  a  very  extensive 

literature  on  the  differential  capacitance  of  solutions  near  either 

solid  (polycrystalline  or  single  crystal)  or  liquid  (mercury) 

electrodes  which  we  will  not  try  to  cover^.  We  should  mention  the 

recent  work  on  the  influence  of  the  crystallographic  orientation  of 

silver  on  the  potential  of  zero  charge  of  the  electrodes,  in  which  a 

detailed  mapping  of  the  influence  of  the  crystal  face  on  the 

21 

differential  capacitance  of  the  inner  layer  is  made 

The  complexity  of  the  system  described  by  the  experimental 
methods  defies  any  simple  theoretical  interpretation.  Yet  these  are 

needed  for  the  understanding  of  what  is  actually  going  on  at  the 
charged  interface.  It  is  clear  that  there  are  two  kinds  of  forces  in 

these  systems:  the  long  ranged  Coulomb  forces  and  the  short  ranged 
forces  that  are  at  the  origin  of  the  chemical  bonds  and  are  also 

responsible  of  the  repulsion  between  atomic  cores.  There  are  important 
quantum  effects  at  the  interface  due  essentially  to  the  quantum  nature 
of  the  electrons  in  a  metal.  For  this  reason  we  have  organized  the 
theoretical  discussion  of  the  chapter  starting  with  very  simple  model 

systems  about  which  a  lot  is  known,  and  to  systems  which  are  much 
more  realistic  but  difficult  to  handle  theoretically.  The  emphasis  of 
the  theoretical  treatment  will  be  on  the  structure  functions,  or 

distribution  functions  p(l),  p  (1,2)  ,...  which  give  the  probability 

•  ij 

of  finding  an  ion(s)  or  solvent  molecule(s)  at  specified  position(s) 
near  the  interface.  The  properties  of  the  interface  can  be  calculated 
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from  these  distribution  functions. 

One  of  the  very  interesting  theoretical  developments  of  recent 

years  has  been  the  exactly  solvable  model  developed  by  Jancovici, 

22 

Cornu  and  co-workers  .  This  is  a  two  dimensional  model  at  a  particular 

value  of  the  reduced  temperature,  and  is  particularly  useful  to 

elucidate  the  subtle  properties  of  the  long  ranged  Coulomb  forces.  For 

the  non  primitive  model  with  solvent  molecules  there  is  a  one 

13 

dimensional  exactly  solvable  model"  .  Exactly  solvable  models  serve  as 
benchmarks  for  approximate  theories  and  to  test  exact  and  general  sum 
rules. 

In  this  chapter  we  will  not  discuss  the  solvent  structure, 

because  this  is  a  subject  under  development,  in  our  discussion  we  will 
restrict  ourselves  to  models  in  which  the  solvent  is  a  continuum  of 

dielectric  constant  e.  The  focus  of  this  article  is  the  structure  of 

the  inner  layer  of  the  interphase,  the  layer  of  ions  that  is  directly 
adsorbed  onto  the  metal  surface.  We  discuss  the  sticky  site  model 
(SSM),  in  which  the  adsorption  sites  are  sticky  points  at  the 

interface.  This  model  requires  the  distribution  functions  of  the 

undecorated,  smooth  surface  as  an  input.  The  theory  of  these 
distribution  functions  is  reviewed  in  section  I.  In  section  II  we 

review  exact  sum  rules  for  the  interface.  In  section  III  we  discuss 

the  structured  interface  in  the  SSM. 

I-THEORIES  FOR  THE  SINGLET  AND  PAIR  DISTRIBUTION  FUNCTIONS 

In  recent  years  there  has  been  significant  progress  in  the 

statistical  mechanics  of  inhomogeneous  charged  systems  such  as  the 

metal  electrolyte  interface.  The  real  interface  is  much  too  complex  to 

be  described  by  a  tractable  model,  so  that  simple  models  that  focus  on 

some  of  the  more  relevant  aspects  of  these  systems  are  used.  Most  of 

the  effort  in  the  past  has  been  directed  to  understanding  models  in 

which  the  metal  side  is  an  ideally  smooth,  charged  surface,  with  or 

without  image  forces.  The  solvent  is  either  a  continuum  of  dielectric 

constant  e,  or  hard  spheres  with  embedded  point  dipoles.  It  is 

possible  to  study  models  with  realistic  solvents  near  charged  or 

neutral  surfaces  in  the  absence  of  electrolytes^,  and  less  realistic 

25 

models  such  as  the  models  with  quadrupoles  ,  or  with  sticky 

26 

octupoles  ,  in  the  presence  of  electrolytes,  but  for  the  time  being 

extensive  calculations  have  not  been  made.  The  primitive  model  of 
charged  hard  spheres  near  a  charged  hard  wall,  all  embedded  in  a 
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continuum  dielectric  is  by  far  the  most  studied,  and  to  date,  the  best 
understood  of  the  three  dimensional  models  of  the  electric  double 

layer. 

There  has  been  considerable  progress  in  the  theory  of  the 

primitive  model  of  the  electric  double  layer.  The  theories  are 

constructed  as  direct  improvements  on  the  original  Gouy-Chapman 

theory,  such  as  the  Modified  Poisson  Boltzmann  (MPB)“  ,  theories,  the 

28 

theories  based  on  functional  density  expansions  and  then  there  are 

the  theories  derived  from  integral  equations.  The  MPB  theories  have 

29 

been  adequately  reviewed  in  another  review  in  this  series  The 

integral  equation  theories  can  be  classified  into  three  groups:  those 

derived  from  the  Omstein  Zemike  (OZ)^  equation  and  those  derived 

from  the  Bom-Green- Yvon  (BGY)^,  Wertheim-Lovett-Mou-Buff  (WLMB)^  and 
33 

Kirkwood’s  equations. 

As  it  happens  in  other  areas  of  condensed  matter  physics, 

the  forces  between  molecules  is  not  known,  but  rather  the  properties 
that  are  measured  are  interpreted  in  terms  of  these  forces.  Since  the 
equations  are  only  approximate,  it  is  impossible  to  estimate  from  the 
comparison  of  the  theory  to  experiment  the  accuracy  of  the  theories. 

For  that  reason  computer  experiments  play  an  important  role,  since 
here  one  knows  the  input  intermolecular  potential,  and  the  properties 
and  structure  can  be  computed  for  systems  consisting  of  a  few  hundred 
to  a  few  thousend  particles.  There  are  two  techniques  used  to  perform 
these  simulations:  molecular  dynamics  (MD),  in  which  the  equations  of 
motion  of  the  molecules  are  solved  simultaneously,  and  the  Monte  Carlo 

method  (MC),  in  which  an  equilibrium  ensemble  is  generated  from  a 

random  walk  algorithm  (a  Markov  Chain).  For  systems  of  hard  charged 
spheres  the  MD  technique  cannot  be  used  because  the  forces  are  both 
singular  and  very  long  ranged,  and  this  produces  insurmountable 
technical  problems.  The  MC  techniques  are  more  amenable  to  simulate 
the  primitive  model  of  the  electric  double  layer^,  although  because  of 
the  long  range  of  the  Coulomb  forces  the  relatively  small  size  of  the 
system  poses  problems'^. 

Models  in  which  the  solvent  is  represented  by  dipolar  hard 
spheres  or  even  higher  multipolar  solvents  have  not  been  simulated  in 
the  neighborhood  of  charged  walls.  The  full  solution  of  the 
HNC1  equation  has  not  been  achieved  for  the  planar  wall,  only  in  the 
the  Mean  Spherical  Approximation  (MSA)  and  in  the  Generalized  Mean 
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Spherical  Approximation  (GMSA)  have  been  solved. 


BASIC  DEFINITIONS:  THE  GOUY  CHAPMAN  THEORY 

We  have  a  mixture  of  ions  of  density  p,  charge  e=ez, 
where  e  is  the  elementary  charge,  z  is  the  electrovalence,  the  ion 
diameter  is  o,  the  number  density  profile  of  i  at  a  distance  z  is 
p(z)=p(l)  from  the  electrode,  which  is  always  assumed  to  be  flat  and 
perfectly  smooth,  (see  figure  1). 

The  singlet  distribution  function  is 
g  (z)=h  (z)  +1 
=p(z)/p , 

where  p  is  the  bulk  density  of  species  i.  The  charge  density  q(z)  is 
given  by 

q(z)=£  ep(z)  (2) 

l 

where  m  is  the  number  of  ionic  species.  The  electrostatic  potential 
pfz)  is  obtained  by  integration  of  Poisson’s  equation 


V20(z)  =  -4jt  q(z)/e 


we  have 


<|>(z)=47t  dt(t-z)q(t)/e 


The  total  potential  drop  Ap  is  obtained  from  (4)  by  either 
letting  z=0  or  z ->-<»,  depending  on  the  reference  potential  of  the 

model.  In  general  the  latter  choice  is  adopted.  An  important  quantity 
is  the  differential  capacitance  Crf  which  defined  by 

Cd=dq5/dA<)  (5) 

where  q?  is  the  surface  charge  on  the  electrode.  This  quantity  is 

seldom  measured  directly,  it  is  inferred  from  either  surface  tension 
measurements,  or  frequency  dependent  AC  measurements  of  the  capacitance. 

The  surface  charge  satisfies  the  electroneutrality  condition 
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oo 


(6) 


dz  q(z) 
o 


=  E  (0)£/4?r 

l 


where  E  (0)  is  the  external  or  applied  field. 

1  sa 

For  low  applied  field,  the  Gouy  -Chapman  -Stem  theory  (or 

modified  Gouy-Chapman  Theory,  MGC)  is  accurate  for  the  primitive  model 

in  comparison  to  computer  experiments. 

Consider  the  Poisson  equation  (3).  If  we  approximate  the  density 

of  the  ions  by  Boltzmann’  distribution  formula 


p(l)=  pexp[-(3e0(l)l 

replacing  into  (3)  we  obtain  the  Poisson  Boltzmann  equation 


(7) 


V20(1)=-4tc/£  £  e(  pexp[-(3e0(l)]  ^ 

1 

This  equation  has  been  solved  analytically  for  several  cases  of 
interest.  A  first  integral  can  be  obtained  multiplying  by  V0(1),  both 
sides  of  this  equation.  For  planar  electrode  this  yields 


E]  =[V0(l)]2=-4trkT/£  Y  p  {exp[-pe<t>(l)]-l}  (9) 

£  Mil 

1 

Using  the  definition  of  Crf  and  the  electroneutrality  relation  (6)  we 
get  the  formula  for  the  differential  capacitance 


C 

d 

where  0(O)=0(z)  |  ^=£0  is  the  potential  at  the  origin,  which  is 
equivalent  to  the  total  polarization  potential  of  the  electrode. 

At  this  point  it  is  convenient  to  make  a  variable  change 

X=exp[-(3e  0(0)1  (n 

where  e  =  ze,  and  z ■  is  the  electrovalence  of  species  i. 

We  integrate  equation  (9)  to  get 


& 


Yep  exp[-(Je  0(0)] 

Mil  1 


(exp[-(3e  0(0)]- 1 } 


(10) 
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where  for  the  general  case,  xo  and  A  are  integration  constants  . 

Since  the  electrovalences  Z(  are  always  small  numbers,  the 

integration  of  the  left  hand  side  is  always  possible  in  terms  of 
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elliptic  functions  For  the  z  =-z,=l  case  the  radicand  of  the  left 
hand  side  is  a  perfect  square  and  the  integral  can  be  performed 
explicitly. 

The  result  for  this  case  is,  for  the  potential  drop  A<t> 

f3eEo/KE=2sinh[A0e(3/2j  (13) 

The  density  profile  is  given  by 


p(z)=p{[l+zae'Kzl/[l  z  a  e_lcz] 

11  l  I 


with  a=tanh[A0  e|3/4].  k  is  the  Debye  parameter 


iC=47t/kTe 


Ip.e“ 


There  are  several  remarks  about  the  Gouy-Chapman  theory. 
From  equations  (10)  and  (14),  we  see  that  for  any  mixture  of  ions  of 
equal  size  the  contact  theorem  is  satisfied. 


kT  £  pX0)=  e/8rc  E*(0)+kT  £  p. 


Contrarily  to  what  has  been  assumed  in  the  literature,  the 

contact  theorem  is  not  satisfied  for  mixtures  of  unequal  size  ions. 
This  can  be  seen  from  equation  (16):  In  the  contact  theorem  all  the 
ions  have  to  be  in  contact  with  the  wall  simultaneously.  Since  there 
is  only  one  distance  in  (16),  there  is  no  way  that  the  theorem  will  be 
satisfied.  However  for  very  high  positive  or  negative  fields  when  only 
one  kind  of  ions  is  present,  and  then  the  Gouy-Chapman  theory  will 
satisfy  asymptotically  the  all  important  contact  theorem. 
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The  density  profiles  obtained  from  the  Gouy-Chapman  theory 
are  monotonous,  that  is  they  show  no  oscillations.  When  the  contact 
theorem  is  satisfied,  and  also  the  electroneutrality  integral  comes 
out  to  be  correct,  p(l)  is  pinned  at  the  origin,  and  has  a  fixed 

integral,  so  that  the  density  profile  cannot  deviate  too  much  from  the 
correct  result.  However,  when  the  density  is  high  the  profiles  will  be 

oscillatory,  and  we  should  expect  deviations  from  the  GC  theory.  This 
is  also  true  for  the  non-primitive  model  in  which  the  solvent  is  a 
fluid  of  finite  size  molecules. 

The  Gouy  Chapman  theory  has  been  solved  for  non-equal  size^  ions. 
For  l-l  electrolytes  the  comparison  to  computer  simulations  shows 
good  agreement,  because  in  the  regime  of  high  electrode  charge  there 
is  only  one  ion  present  in  the  double  layer.  The  situation  should  be 
different  when  we  are  dealing  with  a  mixture  of  anions  or  cations  of 
different  sizes. 

INTEGRAL  EQU ATION S -THE  PRIMITIVE  MODEL 

The  classic  Gouy  and  Chapman^0  theory,  is  based  on  the 

Poisson  equation  and  a  closure  given  by  Boltzmann's  equation.  In  this 
theory  the  ions  in  the  double  layer  are  point  charges  with  no 

exclusion  volume,  all  embedded  in  a  continuum  dielectric.  There  have 
been  a  large  number  of  papers  dealing  with  ways  to  improve  this 
equation.  We  must  remark,  however  that  at  least  in  the  regime  of  low 
density  and  high  temperature  (or  large  dielectric  constant)  the  GC 
theory  is  not  really  bad  in  spite  of  its  simplifications  because  it 
does  satisfy  the  contact  theorem  (16)  asymptotically  for  E-*»,  and  of 
course,  it  also  satisfies  the  electroneutrality  condition  (6). 

However,  since  the  electrochemist  is  really  interested  in  systems  in 

which  the  solvent  is  not  a  continuum  the  primitive  model  should  not  be 
considered  a  working  model  for  the  interpretation  of  experimental 
data,  but  rather  a  learning  model  for  the  theoretician  because  of  the 

availability  of  computer  experiments.  Indeed  as  the  density  and 
coupling  constant  increase,  significant  deviations  from  the  behavior 

predicted  by  the  GC  theory  occur. 

All  of  the  theories  can  be  formulated  as  an  (integral)  equation 
for  the  density  profile  p  (1),  or  alternatively  as  a  (differential  or 
integrodifferential)  equation  for  the  potential  0(1)  Consider  first 
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the  one  particle  direct  correlation  function,  which  will  be  the 
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central  quantity  of  our  discussion  ,  from  which  the  integral 
equations  can  be  deduced: 

cOlsinfp^DAjl+PuO)  (17) 

where  c  (1)  is  the  one  particle  direct  correlation  function,  %.  is  the 

fugacity  of  species  i,  and  u  (1)  is  the  external  potential. 

The  function  c(l)  is  a  member  of  the  family  of  direct 
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correlation  functions  c(l,2..)  ,  which  is  the  sum  of  all  irreducible 

graphs  with  density  factors  p(l)  for  every  field  point  (For  a 

detailed  discussion  of  correlation  functions  see  for  example  Hansen 
and  McDonald^). 

ORNSTEIN-ZERNIKE  BASED  APPROXIMATIONS 
At  the  interface  between  an  electrode  and  a  fluid  the 

density  of  the  fluid  is  a  function  of  the  distance  of  the  point  to  the 

surface  p  (z).  The  Omstein  Zemike  equation  for  this  system  can  be 

obtained  from  an  homogeneous  mixture  in  which  there  are  some  large 
ions,  of  radius  R  -»  «>,such  that  p  R  3-»  0  .  In  this  limit  the  planar 

w  w  w  1 

HAB  (Henderson- Abraham-Barker)  OZ  equation  is4' 


Ml)  -c*(l)  =  P,  [d2  c®<lr,2l)  h,<2>  <18> 

j  ■* 

where  h  (1)  is  defined  by 

I 

h(l)=  [g.(l)-I]-[p.(l)-p.]/p.  (19) 

The  function  h(l)  is  the  density  profile  function  for  ion  i.  The 
function  cw(l)  is  a  much  more  complicated  object,  and  in  general  does 
not  admit  a  simple  diagram  expansion.  To  get  some  insight  about  the 
meaning  of  this  function  we  use  functional  series  expansion^®. 

Consider  the  functional  power  series  expansion  of  ln[p(l)J 
around  the  uniform  density  p. 


Pu.  ( I  )+ln[p .  ( 1  )]=ln[p .  ] + 


+  1/nlV  p  p  ..  d2d3..cB  (1,2,3  ..  )  h  (2)  h(2). 

/  rj  rk  1  jk  J  k 


jk  .  . 


the  direct  correlation  functions  are  defined  by  the  functional 


(20) 
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derivative 


<^(1,2,3...  )=  5nci(l)/6pj(2)8pk(3).. 

The  superscript  B  stands  for  the  bulk  functions.  If  we 
introduce  the  new  function  c*(l),  defined  by 

c*(l )=  -3u  (1 ) - ln[g  ( 1 ))  +h i ( 1 ) 


+  I/n 


!Z p,  pk-f' 


d2d3..c  .  (I ,2,3  ..)  h (2)  h  (3) 

i  jk  J  k 


where  the  inhomogeneous  potential  is  of  the  form 

u(l)=  u,f(l)+w(l) 

I  I  1 

where  u"(l)  is  short  ranged  and  for  a  hard,  smooth  charged 
The  electrostatic  part  is 

w(l)=-e E  z/2 

i  i  i  i 

Combining  this  definition  with  the  functional  expansion 
get  the  HNC1  equation  for  the  flat  wall 


electrode. 


-Pw.(l)-ln[g.(l))=^p.  d2c®  ( |ri2|  )h(2) 
« 


Equation  (25)  has  a  deceivingly  simple  aspect,  but 
because  of  the  long  range  of  w(l)  is  not  convergent,  and 
not  amenable  to  numerical  solution.  Using  the  equation 

c“(lrul)  =  •H(|ru|)  +  c"(lri*>) 


(20)  we 
electrode 


in  fact 
therefore 


w..( )  r  | )  =  ee/er 

ij  1  U 1  i  J  12 

and  replacing  into  (25)  yields 


W>.(l)+ln[g.(l)l-  Y_  P,  W  «”(|ru|)  h.<2> 


where  0(1)  is  defined  by 


0(1)=  E  z 

1  l 


.♦  Y. e.  l62  p> 


(2)/e  r, 


n 


This  equation  is  now  completely  defined  in  terms  of  short 

ranged  quantities,  which  is  not  the  case  for  the  first  form  of  the 
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equation  (25)  .  One  important  observation  about  the  HNC1  is  that  it 
does  not  satisfy  the  contact  theorem  (16),  but  rather^® 

kT  V  p(0)=  e/8re  E2(0)+  p  aP/apA  (30) 

Ld  l  Z  U  0 

with 


P„=  I  P,  <3D 

l 

However,  the  HNC1  satisfies  automatically  the 
electroneutrality  relations  and  the  Stillinger  Lovett  sum  rules,  (see 
below).  The  HNC  is  the  most  accurate  theory  for  bulk  electrolytes.  It 
is  the  theory  that  has  the  closure  with  the  largest  number  of  graphs. 
One  would  expect  that  this  fact  would  remain  true  in  the  plane 
electrode  limit.  However,  because  of  the  inaccuracy  of  the  HNC  for 

uncharged  hard  sphere  fluids  the  HNC1  does  no  do  well  in  representing 
the  exclusion  volume  of  the  ions,  and  is  not  on  the  whole,  such  a  good 
approximation  for  the  electric  double  layer. 

For  high  fields  and  low  concentrations  the  fact  that  we  get 
the  compressibility  rather  than  the  pressure  is  not  very  important  and 
the  HNC1  is  still  a  reasonably  good  theory,  as  will  be  shown  below. 
However  for  dense  systems  this  is  a  rather  severe  shortcoming. 
Specifically,  when  we  are  dealing  with  a  molecular  (dipolar)  solvent 
the  density  is  very  large  and  the  dielectric  constant  e  is  of  the 

order  of  one  (instead  of  80  in  water)  which  makes  the  electrostatic 
term  in  (16)  small  in  comparison  to  the  contact  density  term.  The 

consequence  is  that  the  HNC1  will  put  more  counter  ions  near  the 

electrode  than  the  exclusion  of  the  hard  cores  will  permit.  Eventually 

thermodynamic  stability  conditions  will  be  violated,  such  as  the 

Bogoliubov  inequality^.  For  very  high  charges  the  HNC1  predicts  a 
decreasing  potential  drop  <j)  with  increasing  applied  external  field 
E  52 

Z 

The  bulk  direct  correlation  function 

c®<|r,2l>  (32) 

which  should  be  used  in  solving  the  HNC1  equation  (28)  is  that 

obtained  from  the  solution  of  the  bulk  HNC  equation  for  the  same 

system.  This  however  yields  poor  results  when  compared  to  computer 
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simulations.  Generally  better  results  are  obtained  if  instead  of  the 

HNC  bulk  direct  correlation  function  the  corresponding  MSA  functions 

are  used,  the  general  agreement  with  computer  simulations  improves^. 

There  have  been  several  calculations  with  improved  versions  of 

this  HNC1  equation.  Including  the  next  term  in  the  expansion  (22) 

amounts  to  including  the  bridge  diagram  in  the  bulk  HNC1  calculation. 
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This  was  done  by  Ballone,  Pastore  and  Tosi,  with  good  success  .  The 
density  profile  for  the  1M,  1-1  electrolyte  at  a  surface  charge 

o  =qs<r7e=  0.7,  which  will  be  the  test  case  used  for  comparisons.  This 
is  the  highest  surface  density  simulated,  and  shows  charge 
oscillations  due  to  the  hard  core  of  the  electrolyte.  In  this 
calculation  the  bridge  diagrams  were  computed  directly  from  their 
definition. 


c*(l)=  -(3u  (D-lnfgd)]  +h(l) 


1/ 


pk 


i  jk 


d2d3  h?j(12)  h?k(13)  h®  (32)  h(2)  hfc(3) 


(33) 


where  the  product  of  the  three  bulk  pair  correlation  functions  in  the 
second  term  of  the  right  hand  side  is  first  term  in  the  density 
expansion  of 


cB  (1,2,3)  (34) 

i  jk 

The  results  of  this  calculation  are  shown  in  figure  2.  Since 

there  are  no  adjustable  parameters,  the  agreement  is  very  good.  An 

alternative  less  laborious  procedure  was  suggested  by  Rosenfeld  and 
Blum^,  but  actual  calculations  were  not  performed. 

Another  way  of  inmproving  the  HNC1  approximation  was 
introduced  by  Forstmann  and  co-workers^.  [n  their  method  the  HNC1 
equation  is  used  as  described  above,  but  instead  of  taking  the  bulk 

direct  correlation  function,  as  prescribed  by  (28),  a  local  density 

D 

dependent  c  (r,p)  is  taken.  The  local  density  is  defined  by 

x+A  z +o/2 

P,<x,=  Tar!  dz  I  dz  p.w  <35) 

x  -  A  z  -  o/2 


where  o  is  the  diameter  of  the  ion  and  A  is  an  adjustable  parameter. 
The  bulk  correlation  function  is  then 
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*  * 


c:/ir12i>  p=p  <36> 

■  '  * 

For  the  test  case  with  surface  charge  o  =.7,  the  results  of  this 
method,  as  shown  in  figure  3  are  very  good. 

BGY  BASED  APPROXIMATIONS 
The  BGY  equation  can  be  derived  from  the  one  particle  direct 

correlation  function  c(l).  Consider  again  (20):  letting  the  gradient 
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V  act  on  the  f-bonds  of  the  graphical  expansion  of  c  (1)  »  we  get  the 

BGY  equation 

-kT  V  p(l)=p(l)V  U(l)+y1d2p  (1,2)V  u  (1,2)  (37) 

>•  >  i  >  >j  i  u 

Using  (23)  and  (26)  to  eliminate  the  long  ranged  terms,  we 
obtain  (37)  in  a  different  form 

-kTV  o(l)=p  (1)V  u5r(l)+  p.(t)e  0  (l) 

I  I  l  It  1  1  ll 

r-f  ->  _C  (38) 

d2pt  02 )  V  (  u’ ' ( 12 )  +p .  ( 1  )e  M d2pf2)e  h ..( 12)Vt(  1/er  1 2) 

j  J  j 

This  equation  can  be  integrated  from  °o  to  z,  to  yield 

ln[gi(z)l=-eil0(z)+v(z))+Ji(z)  (39) 

which  together  with  Poisson  equation  (3)  forms  a  closed  system  of 

equations  that  is  very  convenient  for  numerical  solutions. 

This  equation  is  of  the  same  type  as  the  one  derived  from 

the  HNC1  equation  (28):  The  right  hand  side  term  consists  of  three 

contributions:  The  potential  0(1),  which  is  determined  by  the  single 

particle  distribution  function  p.(I),  and  the  terms  Vf{z)  and  J(z) 

which  are  functions  of  the  pair  distribution  function  h  (1,2).  From 

u 

(38)  we  get 

°o 

VXz)  =|dZi  Pi(l)  ^|d2pj(2)ejh  Tl^VG/er^) 

r  (4< 

JU)  =  dz,  Y\d2  P,/12)7,0”0’2) 


We  remark  that  in  (38)  (and  also  in  (39)),  if  the  fluctuation 
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terms  J(z)  and  i|/(z)  are  neglected,  then  we  get  back  the  Gouy  Chapman 

i  ■ 

theory,  which  has  a  known  analytical  solution  . 

In  the  BGY  based  theories  the  pair  correlation  function 
h  (1,2)  must  be  given  by  some  approximation.  The  interesting  feature 
of  the  BGY  equation  is  that  for  no  matter  which  closure,  the  contact 
theorem  (16)  is  satisfied.  However,  the  electroneutrality  conditions 
(see  below)  are  not  satisfied,  in  general.  The  simplest  approximation 
one  could  think  which  is  equivalent  to  Kirkwoods  superposition 
consists  in  writing 

h  (l,2)=hB(r  )  (41) 

ij  ij  i: 


B 

where  h  (r  )  are  the  bulk  pair  correlation  functions,  fails  to 

i  j  12 

satisfy  the  important  electroneutrality  condition  (see  below), 

-e=^|d2pj(2)ejhj(l,2)  (42) 

and  gives  very  poor  results  when  compared  to  the  computer  simulations. 
The  approximation^  ^ 


where 


h  (12)  =  f  (l)f  (2)hB  (r) 
■j  i  j  i  j 


r  >C 

12  ij 


=-l 


r„<0 

12  ij 


a  =(a+o)/2 

U  *  j 


(43) 

(44) 


is  constructed  so  that  the  functions  f.(z)  are  required  to  satisfy  the 
electroneutrality  condition  (42)  for  the  inhomogeneous  pair 

distribution  function.  This  yields  an  integral  equation  for  those 

functions.  This  construct  of  the  inhomogeneous  pair  correlation 
function  can  give  negative  values  of  gTl.2).  The  problem  is 
specially  severe  for  dense  systems,  beyond  2  M  of  salt  concentration 
or  also  for  the  test  case  of  high  surface  density  charge.  A  simple  way 
to  circumvent  this  problem  was  suggested  by  Caccamo,  Pizzimenti  and 
Blum59 
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(45) 


h  (1.2)=  f  (l)f (2)hB (r)+A  (1.2) 

i)  1  j  i  j  >j 


r,2>0 


=-l  r  <o 

i-  >j 

where 

A  (l,2)=AhPY(r  )  z.z  <2o  (46) 

i)  i  j  12  I  2  ij 

the  parameter  A  is  adjusted  to  eliminate  negative  correlations  from 
h  (1.2).  The  results  are,  however  not  very  sensitive  to  the  exact 

'j 

value  of  A.  The  comparison  to  the  Monte  Carlo  simulations  is  again 

*  -i 

good,  this  method  yields  for  the  test  case  with  a  =EJeo~=J  the 

density  oscillations  in  the  profile  of  the  counterions.  Figure  4  shows 

34 

the  comparison  to  the  computer  simulations  of  Torrie  and  Valleau 
The  inhomogeneous  pair  correlation  function  h  (1,2)  can 
be  obtained  from  the  inhomogenous  OZ  equation 


h  (l,2)-c  (1,2) 
u  u 


I 


d3  hik(l,3)pk(3)  ckj(3,2) 


and  a  suitable  closure  for  the  direct  correlation  function. 


(47) 


c  .(l,2)=-(3  u  (1,2)  rr>0,j  (MSA2  approximation)  (48) 

Cj(l,2)=-P  u.(l,2)  +  htj(l,2)-ln(g|j(  1,2)1  (49) 

r  >a 

12  i} 

(HNC2  approximation) 

The  MSA2  approximation  cannot  be  integrated  explicitly,  as 
is  the  case  of  the  homogeneous  MSA.  When  the  ions  are  approximated  by 
charged  points  then  for  some  specific  form  of  the  density  profiles 
p(l)  the  OZ  equation  can  be  integrated,  and  series  solutions  have 
been  given  (Blum,  Hernando  and  Lebowitz^  and  Camie  and  Chan^).  The 
numerical  solution  of  the  MSA2  and  HNC2  has  been  extensively  studied 
by  Henderson  and  Plischke*^.  Approximate  solutions  of  this  kind  will 
satisfy  the  contact  theorem  (16)  as  well  as  the  electroneutrality  and 
dipole  sum  rules*^ 

The  BGY-HNC2  equation  has  been  solved  numerically  by 
Nieminen,  Ashcroft  and  collaborators*^,  however  for  systems  with 
neutral  molecules. 
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WLMB  BASED  EQUATIONS 

Yet  another  integral  equation  is  derived  from  the  one 
particle  direct  correlation  function  c  ( 1 )  (22)  by  introducing  relative 
coordinates  for  all  field  in  the  diagram  representation  and  taking  the 
derivatives  with  respect  to  those  coordinates^  This  yields  and 
exact  hierarchy  of  equations  that  is  related  to  thy  BGY  hierarchy.  The 
first  member  of  the  Wertheim-Lovett-Mou-Buff  (WLMB)  equation  is 

Vp(l)+  (3  p(l)  Vu(l)  =  pd) 

This  equation  contains  long  range,  divergent  terms.  Introducing  the 
local  potential  0(1)  (29),  we  have^  ^ 

Vln[p(l)]+p  V<|>(1)+  p  Vu"(l)  = 

66 

This  equation  has  been  studied  by  Henderson  and  Plischke  in 

detail.  As  can  be  seen  in  the  figure  5,  it  yields  very  good  results 

* 

for  the  test  case  of  a  =.7.  The  calculations  were  performed  solving 
both  the  HNC2  closure  for  the  inhomogeneous  pair  correlation  function, 
and  also  the  MSA2  closure  in  a  few  cases. 

A  simplified  version  of  the  WLMB  equation  that  produces 
reasonably  good  results  was  studied  by  Colmenares  and  Olivares^ 

KIRKWOODS  EQUATION 

An  interesting  approach  has  been  suggested  by  Kjellander  and 
Marcelja  ,  based  on  the  observation  that  for  the  HNC  approximation  the 
chemical  potential  can  be  obtained  explicitly  as  a  function  of  the 
pair  potential  h..(ri2)  for  an  homogeneous  fluid.  Then,  within  the  HNC 
the  function  c.(l)  can  be  explicitly  evaluated. 

The  central  idea  is  to  slice  the  three  dimensional  space 
into  two  dimensional  layers  that  are  homogeneous.  The  three 

dimensional  OZ  equation  can  be  mapped  into  coupled  set  of  N  two 
dimensional  OZ  equations  for  a  mixture  of  N  components,  each  component 

is  an  ion  in  a  different  layer.  The  particles  interact  with  a  species 

dependent  interaction  pair  potential.  In  the  limit  of  an  infinite 

number  of  layers  this  procedure  yields  the  correct  inhomogeneous  OZ 
equation.  The  chemical  potential  p(a)  of  the  ith  ion  in  the  ath 


3 


d2  c  sr(  1,2)  V  o  (2) 

>J  2  J 


(51) 


d2  c  (1.2)  V  p(2) 


(50) 
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layer  is  given  by  Kirkwoods  equation: 

^(a)  =kT  In  P((a)+  kT  In  Aq/Az+  V(  (a) 

fi  -  (52) 

+  V  p  (13)  dX  dR  g  (R,aP;X)  9  g,j ( R,aP;A,) 

(3,7  0  ** 

where  X  is  the  coupling  parameter,  Az  is  the  thickness  of  the  layer, 

Aq  is  the  ideal  gas  fugacity,  V(a)  is  the  interaction  between  a 

panicle  in  layer  a  and  the  wall;  R  is  the  two  dimensional  distance. 

In  the  HNC  closure 

c  (R,a{3)=  h  (R,a(3)-u  (R,a(3)AT-ln[g  (R,aP)J  (53) 

ij  y  ij  ij 

Kirkwoods  equation  can  be  integrated  to  yield 

r 

p(a)  =Az/AQexp-  Pp(a)  +  Vp  ( (p)  dR["l/2h2  (R,aP)-ci  (R,ap) 

l  '  PT7  '  J 

-Uj(R,ap)/kT  -1/2  In  ^.(R.aPi/explu.^R.aP) /kT]JR=0  -  O.(a) 

(54) 

where  O(a)  is  the  average  potential  for  layer  a. 

<t>.(a)  =  2xe2/e  £  P,(P)|za-zpl  (55) 

P 

★ 

The  results  for  the  same  case  of  o  =  7  is  shown  in  figure  6.  The 
agreement  is  also  excellent. 
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II-  SUM  RULES  FOR  THE  CHARGED  INTERFACE 
The  sum  rules  for  the  charged  interfaces  can  be  classified  in  two 
categories: 

1)  the  dynamic  sum  rules  which  are  derived  from  balance  of  forces 
considerations, 

2)  the  screening  sum  rules,  which  are  specific  to  Coulomb  forces. 
Because  of  the  very  long  range  of  the  electrostatic  forces,  the 
stability  of  the  system  requires  that  all  charges  surround  themselves 
with  a  neutralizing  cloud. 

DYNAMIC  SUM  RULES 

Consider  the  Bom  Green  Yvon  (BGY)  equation"^  ^ 


-kT  ViP.(l)«  p(l)  Vu(l)  +)  ]d2  ptj(l,2) 


V  u  (1,2) 
i  ij 


(56) 


where  k  is  Boltzmann’s  constant,  T  is  the  absolute  temperature,  pX  1 ) 
is  the  density  of  species  i  at  position  r  sx  .y^Zj.  If  species  i  is 
non  spherically  symmetric  then  the  integration  should  also  include  the 
rotational  coordinates,  the  Euler  angles  u  (1)  is  the 

single  body  potential  acting  on  i.  The  pair  potential  is  represented  by 
u  (1,2).  Similarly  p..(l,2)  is  the  pair  density  function,  which  can  be 
written  in  the  form 

p  (l,2)=p(2)p(2)g..(l,2)  (57) 

where  g  (1,2)  is  the  pair  distribution  function, 
y 

For  charged  particles  it  is  convenient  to  separate  the 

electrostatic  from  the  nonelectrostatic  contributions.  For  the  singlet 

potential  we  have 

u(l)=e<t>(l)+u”(l)  (58) 

1  I  t 

where  <tK  1)  is  the  local  potential  at  r ,  e  is  the  electric  charge  of 

species  i  and  u”(l)  is  the  sum  of  all  nonelectrostatic  forces,  which 

includes  the  hard  repulsion  from  the  wall.  Similarly,  the  two  body 

contribution  is 

V1-2)=','/eVu"(l'2)  (59) 

where  u"(l,2)  is  the  nonelectrostatic  interaction  of  the  pair  ij,  and 

•j 

the  first  term  is  the  electric  potential  between  the  charges  of  i  and 
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j.  Furthermore  we  use  the  function 


and  the  potential 


h  (l,2)=g  (1,2)-1 

ij  ®ij 


V  (l)=<)>(l)-t-^e|ci2pa)/eri2 


where  r  =|r  -rj.  Thc  BGY  equation  is  then  transformed  to 


-kTV  p  (l)=p  (1)V  usr(  l)+p  (l)e  V  w  (1) 

ir i  r  i  1  i  r  i  iii 

+V|d2p  ( 1 ,2)5  usr(  l,2)+p  ( 1  )e .  V [d2p  (2)e  h .( 1,2)V  ( 1/er  ) 

/  lJ  1  ‘J  >  J  J  iJ  1  12 

J  J  J 

(62) 

We  consider  a  system  which  is  limited  by  an  arbitrarily  rough 
planar,  but  charged  surface.  The  precise  mathematical  requirement  is 
that  there  is  a  prism  with  an  arbitrarily  large  cross  section  area  S  , 
and  height  L  (  the  volume  V=SL),  such  that  the  force  through  the  walls 

Iff 

parallel  to  z  is  of  0(S  °),  where  8>0,  as  S-*». 

We  integrate  now  over  the  whole  prism  of  volume  V.  This  is  a 
generalization  of  the  case  of  a  decorated  surface  recently  discussed  in 

71 

the  literature  ,  and  is  of  interest  in  the  discussion  of  the 
electrochemistry  of  rough  surfaces. 

We  integrate  the  BGY  equation  in  the  volume  of  a  prism  of  the  same 
section  S  but  smaller  height  Lj<L.  Summing  over  all  species  i,  we  get 

-kT  ^|dlVipi(l)  =  ^|dl[p.(l)Viu;r(l)+  p.(l)e.  V^.G )  ] 

i  i  J 

♦£]«  n<*2[  P./1’2>V"(l’2)+'\<l)<V2)  Ml, 2)  Vleeyer^)] 


l  J 


Since  in  our  system  V)=a/az]  the  integral  of  the  left  hand  side 
term  can  be  easily  performed 
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drV^d)  = 


r  r  1 

)  =  dx(dy(  dz(  3p(l)/azt 
s 

"  ^ 

=  S  -  1/S  dx)dyi  p^x^.z) 


driVip[(l)=  S  ptt^)  -  p(S) 


where  we  have  defined  the  average  contact  value  as 

pt(0)  =  1/S  dx,dy ,  P/VVV  (65 

To  integrate  the  second  term  in  the  right  hand  side  we  use 
Poisson’s  equation 

(1)=  -4  7t/e[  epO)  =-V.E(l)  (66) 

where  we  have  used  the  electric  field 

E(l)=-Vy  (1)  (67) 

Substitution  into  the  second  term  of  the  right  hand  side  leads  to 


dr^[  e.p.(l)j  [  Ed)]  =  -  e/4it  dr(  (1)  Ez(l) 


=  e/4jr 


[dr,[v,  .  E(l)] 


E(l) 


where  E(l)  is  the  electric  field  at  position  r,  and  Ez(l)  is  the  z 
component.  Integrating  by  parts,  we  get 


r  f  1 

dri[[  CIP,(1)]  [  Ez(1)1=c/8jc  dx jdyj  dz,  dE\\)/dz 

Js  Jzs(x,,y,)Z 

+  e/4 n  dx  dy  dz  [e  (1)  aE  (l)/ax  +  E  (1)  aE  ( 1  )/ay  1 

i  i  i  I  z  x  i  z  y  •  i 
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The  second  term  of  the  right  hand  side  is  zero:  For  a  periodic 
interface  in  the  x  and  y  directions,  if  we  take  S  to  be  the  surface  of 
a  unit  cell,  then  the  terms  like  aE  (l)/axi  will  be  equal  but  with 
opposite  sign  for  neighboring  cells.  For  the  general  random  interface 
we  conjecture  that  this  term  is  finite:  then  in  the  limit  S-*»  the 
contribution  vanishes.  We  have 


1/S  jdr  [[  ep(l)J  [  Ez(l)]  =-e/8t:  <E^>s 
where  the  average  square  field  in  the  z  direction  is 


<E‘>C=  1/S 

l  O 


tlx.rty,  E-i(x|.yi,z  ) 

s 


The  first  term  yields 


1/S 


dl  p  (1)V  u5r(l)=  <p (1)  ausr(l)/az  >~ 

I  1  I  I  1  13 


(70) 


(71) 


(72) 


where 


<PU)  au”(l)/az  >s=  1/S 


dx)dyi  dz(  p  ( 1 )  du'{\)/di{ 

s  *,(v>v 


(73) 


The  last  term  in  the  right  hand  side  of  equation  (63)  is  the  average 
of  the  pair  forces.  We  write 


!p= 

dx.dy, 

dx^ 

S  J 

s 

= 

dx.dy, 

dx2dy2 

■ 

s  J 

S 

L( 
d  z 


z,(Vy,> 

L 


dZ2 

zs(x2,y2) 


ZF*i<u 


>  J 


1 

d  z 


i 


z  (x,y  ) 


L 

r  i 

dz2  + 


L  Jz  (x,y) 


L 

dz. 


i  J 


2) 


JVa.2) 

i  j 


(74) 


F.,0.2)-  [  Pj/t'2)^Iu"(l,2)+P1(l)p32)  h/1.2)  VCee/er,^] 

(75) 


where 


clearly 

F  (1,2)= -F  (2,1)  (76) 

u  j> 

from  where  the  First  double  integral  in  (74)  vanishes.  Since  the 
interactions  are  short  ranged,  the  second  double  integral  yields 
the  virial  contribution  to  the  bulk  pressure  P 


P  =1/S 


dx,d>', 

S 


dx;dy, 
S  ‘ 


,  i 

dz, 

z  (x,y) 


kp,<‘-2> 

L,  >  j 


(77) 


=  2jt/3  )  p  p  dr  r3  fe  eh 

/  »  j  [  *  j  1  j 

o 


(r)/er'+  g  (r)  au  (r)/ar 
>  j  u 


Putting  it  all  together  yields  the  general  contact  theorem  for  a 
surface  that  is  planar  on  the  average,  but  not  necessarily  smooth: 

kT£  p(0)=  P  +  e/8jt<Ez>s  -  £  <p(l)  au^D/az^  (78) 

1  i 

This  theorem  is  a  generalization  of  the  ’-reviously  derived  contact 
theorems  to  the  realistic  case  of  non  smooth  electrode  surfaces.  It 
contains  the  previous  results  as  particular  ca?-"S.  If  the  interface  is 
a  smooth  hard  wall  then  the  surface  averages  become  the  surface  values 
of  the  parameters  and  we  get 

kT[  pX0)=  P  +  e/8n  E*(0)-  J]  <pO)  au”(l)/az  >  (79) 

1  1 
where  the  last  term  is  now 


<p (1)  ausr(l)/az  >  = 

i  i  i 


dz  p  ( 1)  au!r(I)/az 

i  i  i 


(80) 


o 


when  uf  (1)  is  zero  then  we  recover  the  previously  obtained  results 

'2/ 


kT  l  pX0)=  P  +  e/8 k  E^(0) 


(81) 


for  the  primitive  model  with  a  continuum  solvent  of  dielectric  constant 


23 


e,  and 


kljj  p(0)+  [  Pn(0)]=  P  +  l/3Jt  E\ 


for  the  molecular  solvent  case.  Here  the  second  sum  in  the  left  hand 
side  is  over  all  the  neutral  molecules  in  the  system,  while  the  first 
is  over  all  the  ions  present. 

When  the  surface  has  an  array  of  sticky  adsorption  sites,  such  as 
in  the  case  of  the  SSM  model  discussed  below,  then  the  adsorption 
potential  has  the  form 

exp[-j3ua(r)]=  1  +  >-a(R)8(z)  (g3) 


^a(R)-  y^A»a  8(R-njaj-n2a2) 


Here  R=x,  y  is  the  position  at  the  electrode  surface,  and  z  the 
distance  to  the  contact  plane,  which  is  at  a  distance  a/2  from  the 
electrode.  In  (84),  nj,n2  are  natural  numbers,  and  a^^  are  lattice 
vectors  of  the  adsorption  sites  on  the  surface.  The  number 
represents  the  fugacity  of  an  adsorbed  atom  of  species  a.  Define  now  the 
regular  part  of  the  density  function 

Pi(l)=y.(l)exp[-Pu"(l)]  (85) 

Replacing  into  the  general  contact  theorem  (78)  yields 


kT£  pX0)=  P  +  e/8n<E*>s  +  kTX  [oy^U/az^ 


This  theorem  has  been  verified  recently  by  F.Comu,  for  the 
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exactly  solved  model  of  a  one  component  plasma  in  two  dimensions 
THE  SCREENING  SUM  RULES 

In  systems  that  are  electrically  neutral,  any  fixed  arrangement  of 
charges  is  screened  by  the  mobile  charges  of  the  system.  In  homogeneous 
bulk  phase  this  is  an  intuitively  natural  fact,  because  if  the  long 
ranged  Coulomb  forces  would  not  be  screened  then  the  partition  function 
would  not  exist  (it  would  diverge),  and  matter  would  not  be  stable. 
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This  is  expressed  by  the  fact  that  the  charge  distribution  around  a 
given  charge  e|  is  of  equal  value  but  opposite  sign. 

In  the  homogeneous  bulk  phase  this  is  a  natural  fact: 

'c‘=Z}d2pj(2)  *  h‘j(1,2)  (87) 


Rotational  invariance  in  bulk  fluids  require  that  not  only  charges 

but  also  multipole  of  arbitrary  order  should  be  screened  by  the  mobile 
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charges  of  the  media  This  fact  is  much  less  intuitive  in  the 

neighborhood  of  charged  objects,  in  particular  in  the  neighborhood  of  a 

charged  electrode.  However  the  theorems  hold  and  in  classical 

mechanics,  at  least,  perfect  screening  of  all  multipoles  occurs,  in  the 

homogeneous  or  inhomogeneous  systems. 

These  conclusions  are  supported  by  the  results  of  the  exactly 

solved  Jancovici  model7^.  However,  perfect  screening  of  all  multipoles 
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does  not  occur  in  quantum  systems  or  in  systems  out  of  equilibrium 

As  a  consequence  of  the  screening  the  second  moment  of  the  pair 
distribution  function  must  be  normalized.  This  is  the  Stillinger-Lovett 
moment  relation7**.  As  was  shown  by  Outh waite77,  it  can  be  written  in  the 
form  of  a  normalization  condition  for  the  electrostatic  potential 

¥,(r)=  1/e  [e/h£  pJdr,  ^(r, V<  |  r-r,  | )]  (gg) 

J  J 

We  have  then  that  the  second  moment  relation  is 


1/kT 


(89) 
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Camie  and  Chan  have  shown  that  this  normalization  condition  is 

also  valid  for  the  inhomogeneous  systems  of  charged  particles:  Consider 
the  inhomogeneous  Omstein-Zemike  (OZ)  equation 


h  (l,2)-c  (1,2) 


(90) 


where  c  (1,2)  is  the  direct  correlation  function  and  the  singlet 
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density  function  p(l)  satisfies  the  electroneutrality  condition 


Y  e,|dri  p,<i)=-c 


where  qs  is  the  surface  charge  density  and  S  is  the  area  of  the 
interface.  The  pair  distribution  function  satisfies  the  sum  rule 


-e=  ][Jd2p(2)  Cj  ha ,2) 


From  the  diagram  expansion  we  write  the  direct  correlation  function  as 


c  (l,2)=c°(l,2)-3uel(l,2) 
u  u  >j 


where 


ue(l,2)=  eeyer 

ij  i  J  12 


is  the  electrostatic  part  of  the  interaction  and  c°(l,2)  is  the  short 

r  ij 

ranged  part  of  the  direct  correlation  function,  in  the  dense 

media 

0(!.2te.-u;'(r,2>^|dr)u;;(r,3)pk(rJ)hk(3.2)  (95> 

and  from  the  OZ  equation  (90)  we  get 

h..(l,2)=-p*(l,2)  e  -c°(l,2)  +  V|d3  c°  (l,3)p.(3)  h  (3,2)  (96) 

ij  rT|  i  ij  J  ik  k  kj 


Multiplying  this  equation  by  e^pX2),  integrating  over  r2  and 
summing  over  j  ,  yields  after  use  of  the  electroneutrality  condition 


1/kT^eJdr  j  p.(  1  X|>.(  1  t2)=  1 
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which  is  the  generalization  of  the  Outhwaite  formula  of  the  Stillinger 
Lovett  sum  rule  for  inhomogeneous  charged  systems. 

OTHER  SUM  RULES 

For  flat  hard  electrode  surfaces  there  are  number  of  other  sum 
rules.  A  complete  review  of  these  rules  was  recently  made  by 
Ph.A.Mamn^.  We  will  just  mention  a  few  of  the  more  relevant  to  the 

on 

calculation  of  density  profiles  in  the  electric  double  layer 

kT  aln  p.(D/aE0  =  £  d2p.(2)  e.  ha.2)  (z2-z,)  (98) 

j  J 

where  E0  is  the  bare  field  at  the  electrode  surface.  The  differential 
capacity,  which  is  defined  by 


C  =  a  q  /a  0.(D 

d  s  i 


(99) 


where  is  the  surface  charge,  q5=EQe/4rt,  and  0(1)  is  the 
potential  drop,  satisfies  the  sum  rule 


1/Cd  = 


8tc2/e2S  y  jdld2p.(2)  p.(l)  ec.h  (1,2)  (z2-z,)2 

i7  ^ 


(100) 


The  surface  tension  y  obeys  relations  that  can  be  given  in 
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terms  of  the  direct  correlation  function 


y  =(ickT/2)^  dZjdz2  ap.(l)/az)  apj(l)/sz2 


dr,2r;2  c  (U)  (101) 


Using  the  Omstein  Zemike  equation  we  get  the  form  that  contains 
the  pair  correlation  function 


y  =(ji/2kT)^T  dz(dz2  p.(l)au.(iyaz1  p.(2)  au.(l)/az2 


(102) 


These  sum  rules  provide  ways  of  asserting  the  accuracy  of  the 
different  approximations  used  to  compute  the  charge  and  ion  density 
near  charged  walls. 
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Ill-  THE  STICKY  SITE  MODEL 

The  theoretical  discussion  of  smooth  surfaces  is  considerably 
simpler  than  that  of  a  realistic  surface,  in  which  the  solid,  usually  a 
metal,  has  a  well  defined  crystal  structure.  The  reason  is  that  in  the 
case  of  a  smooth  surface  the  problem  is  one  dimensional,  rather  than 
three  dimensional.  The  analysis  of  a  realistic  metal  surface  potential 
in  contact  with  an  ionic  solution  is  extremely  difficult,  and  requires 
the  use  of  very  large  computers.  We  would  like  to  discuss  a  simple 
model  of  a  structured  interface  which  predicts  surface  phase  behavior, 
for  the  adsorbed  layers,  and  which  is  mathematically  tractable.  In 
fact,  if  the  correlation  functions  of  the  smooth  surface  model  are 
known  to  all  orders,  then  the  properties  and  correlation  functions  of 
our  model  can  be  computed  exactly,  at  least  in  principle. 

on  n  j 

The  model8-3  4  combines  two  ideas  that  have  been  used  a  very  long 
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time  ago  :  Boltzmann’s  sticky  potential  ,  and  the  adsorption  site  model 
Rfi  87 

of  Langmuir  .  The  elegant  work  of  Baxter  in  which  the  Percus-Yevick 
approximation  of  the  sticky  hard  sphere  model  is  solved  and  discussed, 
shows  that  this  model  has  a  particularly  simple  mathematical  solution. 
In  Baxter’s  work  the  potential  has  the  form 

exp  [*(3u4(r)]=  l  +  X5(r-o*)  (103) 

where  fi=l/kT  is  the  usual  Boltzmann  thermal  factor,  u(r)  is  the 

intermolecular  potential,  A.  is  the  stickiness  parameter,  r=(x,y,z)  is 
the  relative  position  of  the  center  of  the  molecules,  and  a  is  the 

diameter  of  the  molecules.  The  right  hand  side  term  represents  the 
probability  of  two  molecules  being  stuck  by  the  potential  ua(r):  this 
occurs  only  when  the  two  molecules  actually  touch,  and  for  this  reason 
we  use  the  Dirac  delta  function  S(r-o'),  which  is  zero  when  the 

molecules  do  not  touch,  is  infinity  when  they  do,  but  the  integral  is 
normalized  to  one.  The  stickiness  is  represented  by  the  parameter  X, 
which  except  for  a  normalization  factor,  can  be  considered  as  the 

fugacity  of  the  formation  of  the  pair. 

The  Langmuir  adsorption  sites  can  be  represented  by  a  collection 
of  sticky  sites  of  the  same  form  as  was  suggested  by  Baxter.  Only  that 
now  we  do  not  have  a  sphere  covered  uniformly  by  a  layer  of  glue,  but 
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rather  a  smooth,  hard  surface  with  sticky  points,  which  represent 
adsorption  sites  where  actual  chemical  bonding  takes  place.  For  this 
model,  equation  (103)  has  to  be  changed  to 

exp[-(Jua(r)]=  1  +  Xa(R)8(z)  (104) 

with  — 

*-a(R)=  )  8(R-n1a1-n2a2)  (105) 

nln2 

Here  R=(x,y)  is  the  position  at  the  electrode  surface,  and  z  the 
distance  to  the  contact  plane,  which  is  at  a  distance  a/2  from  the 
electrode.  In  (105),  nj,n2  are  natural  numbers,  and  aj,a^  are  lattice 
vectors  of  the  adsorption  sites  on  the  surface.  The  number  Xa 
represents  the  fugacity  of  an  adsorbed  atom  of  species  a,  onto  the 
surface  which  has  a  perfectly  ordered  array  of  adsorption  sites.  While 
this  is  not  a  requirement  of  the  model,  it  makes  the  mathematical 
discussion  much  simpler. 

Consider  now  a  fluid  consisting  of  only  one  kind  of  particles  of 
diameter  a  ,  near  a  smooth,  hard  wall  with  sticky  sites.  The  fluid  has 
N  particles  and  the  volume  of  the  system  is  V.  The  Hamiltonian  of  the 
system  is 

H=  H°+H2  (106) 

where  H®  is  the  hamiltonian  of  the  system  in  the  absence  of  the  sticky 

2 

sites  on  the  hard  wall,  and  H  is  the  sticky  sites  interaction 

N 

H2  =  )  u,(ri}  (107) 

i  =  l 

where  u^(r)  is  the  sticky  interaction  of  eq.(104). 
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The  canonical  partition  function  of  this  model  is 
1  f  a  N 

Z=  exp  H3HU)  n  n  t  1  +yRi)8(zi)]  dr.  (108) 

»  a-.  i=l  a 
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Expanding  the  product  in  (108)  and  integrating  the  Dirac  delta 
functions  we  get,  using  the  single  component  notation  to  avoid  heavy 
and  unnecessarily  complex  equadons,  with  the  understanding  that  in  the 
multi  component  case  X  is  a  vector  quantity  with  components  X...X,  ,  and 

1  dL 

the  necessary  modifications  of  N!  and  the  integrations  have  to  be  made. 

Z  =Z°^W)p«(rr.  rs,  u09) 

all  sites 
on  the  surface 

where  r=  R^O  is  the  position  of  the  ith  adsorbing  site  on  the 

surface.  p^(i-j...rs)  is  the  s-body  correlation  function  of  the  smooth 
interface  model. 


p“(r,...  rs)  - 


1 


(N-s)!Z 


770 


N 


exp[-(5H]n  dr. 
i=s+l 


=  g?(rr-.rs) 


(110) 

(111) 


Equation  (111)  defines  the  s-body  correlation  function  g®(r....r  ) 

0  sis 

while  p  (r)  is  the  singlet  density  of  the  smooth  wall  inhomogeneous 

problem.  The  smooth  wall  partition  function  is 


70  _1 
z  NT 


jO, 


N 


exp  [-pH  1  nidrj 


(112) 


The  important  observation  is  that  in  our  sticky  sites  model  (SSM), 
the  excess  properties  of  the  interface  depend  only  on  the  correlation 
functions  of  the  smooth  interface  model.  In  fact,  introducing  the 
potentials  of  mean  force  w®(r,...r  ) 

3  1  5 


®s(rl  rs)=  “p[-P  ws(rl-rs)] 


(113) 
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z  /z°.  I  «p[*  E  -x-v] 

s ! 

«1  1  Iltt! 
on  tne  surface 

(114) 

The  left  hand  side  of  (114)  is  closely  related  to  the  excess  free 

energy  of  the  sticky  site  model  as  compared  to  the  smooth  wall 
problem.  In  fact,  from  standard  thermodynamics  the  excess  free  energy 
is 

AFS=  -1/p  ln(Z/Z°)  (115) 

It  is  clear  that  in  (114)  the  right  hand  side  is  the  grand 

canonical  partition  function  of  the  two  dimensional  lattice  gas  in 

which  the  fugacity  is  a  =Xp^(r-)=Xp^(0),  where  we  have  used  the  fact 

that  the  single  particle  density  of  the  smooth  wall  problem  depends 

only  on  the  distance  to  the  electrode  z,  and  not  on  the  position  of  the 

adsorbed  particle  on  the  interface.  The  SSM  is  a  model  that  decouples 

the  structure  of  the  interface  from  the  inhomogeneous  charged  fluid 

problem,  which  in  itself  is  a  very  difficult  one,  even  for  the  smooth 

interface.  The  excess  properties  of  the  lattice  surface  are  formally 

the  same  as  two  the  dimensional  lattice  gas  problems  that  have  been 
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extensively  investigated,  notably  in  connection  with  the  Ising  model 
If  the  area  of  the  interface  is  S  ,  and  the  number  of  sites  is  M,  the 

number  of  sites  per  unit  area  is 

0)  =  M/S  (116) 

Define  also  the  excess  free  energy  per  unit  area 

b?  =  AFs/S  (117) 

and  the  fraction  of  occupied  sites  0  =N  /M  of  particles  of  species  a. 
If  we  compare  (109)  to  the  grand  canonical  partition  function  of  a 
lattice  gas  on  the  sites  of  our  model,  we  see  that  the  number  of  sites 
occupied  by  particles  of  species  a  is 

<Na>=  a  In  Z N/a  In  Xa  (118) 

and  using  (116)  the  fraction  of  occupied  sites  is 

ea=(  l/M)alnZN/aln\a  (119) 
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The  excess  free  energy  Af5  is  also  the  excess  pressure  due  to  the 
presence  of  the  discrete  sites  structure  on  the  interface. 


EXACT  RELATIONS  FOR  THE  SSM 

SUM  RULES  FOR  THE  FLUID  DENSITY  FUNCTIONS 

The  SSM  for  a  charged  interface  satisfies  a  number  of  exact 

relations:  The  first  one  is  the  analog  to  the  Gibbs  absorption 
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equation:  If  we  define  the  surface  excess  in  the  SSM  of  species  a 

oo 

Ar  =f  dz[p  (z)-p^(z)] 

a  J  o  a  a  (120) 

where  p^(z)  is  the  density  profile  of  the  smooth  wall,  and  p  (z)  is 

a  a 

that  of  the  SSM,  then  we  must  have 


Ar  s-flaA^/SlnX 

A  d 


(121) 


Because  of  the  singularity  in  the  f-function,  the  distribution 
functions  also  must  have  a  singularity.  They  must  be  of  the  form 

«  f  1  .  .  r-\  (122) 


pa(r)=[l+yR)5(z)]ya(r) 


where  y(r)  is  the  regular  part  of  the  distribution  function.  The 

2L 

average  number  of  adsorbed  particles  is  given  by  the  integral  of  the 
singular  pan  of  the  density  p„(r) 

a 


<Na>=JVR)8(z)ya(r)dr 


(123) 


since  this  integral  is  actually  a  sum  over  the  sticky  sites  m,  we  get 
the  relation 


W»<rn.> 


(124) 


which  is  exact  :  The  stickiness  parameter  Xa  is  the  ratio  between  the 
number  of  adsorbed  and  non  adsorbed  particles  at  site  rm  of  the 
lattice.  Integrating  (122)  over  all  space  yields 

N=<Na>+Jya(r)dr  (125) 

which  can  be  rewritten  as 

<Na>=-Jdr[ya(r)-pal  (126) 

Because  of  the  analogy  of  the  partition  function  to  the  grand 
canonical  partition  function,  a  relation  similar  to  the  compressibility 
relation  holds 
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(127) 


<N2>-<N  >2=X  a  <Na> 

EXACT  RESULTS  FOR  THE  ADSORBED  LAYER 

Consider  the  case  of  a  simple  salt  dissolved  in  water,  near  a 

metallic  electrode:  in  the  SSM  there  will  be  three  components,  the 

anion,  the  cation  and  the  solvent,  and  the  lattice  atoms.  In  the 

limiting  case  of  the  SSM,  the  sizes  of  the  different  species  play  a 

crucial  role  in  the  possible  ordering  of  the  ad  layers  at  the 

interface.  It  will  be  convenient  to  picture  the  ions  as  having  a  hard 

sphere  core  with  a  diameter  oa,  o^,  and  the  solvent  as  having  a  hard 

core  with  diameter  on-  The  lattice  spacing  of  the  metal  surface  is  d, 

and  because  it  is  the  most  stable  surface,  we  will  restrict  ourselves 

to  discussing  the  (111)  surface  of  the  fee  crystals,  or  the  (100)  face 

of  the  hep  crystal,  that  is,  the  triangular  lattice.  In  the  most 

general  case,  all  three  components  can  be  adsorbed  competitively  and 
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this  situation  can  give  raise  to  very  complex  phase  diagrams  .  This 
most  general  case  can  be  modeled  by  the  spin  S=1  Ising  model,  has  a 
very  rich  phase  diagram,  involving  first  and  second  order  phase 

transitions  and  multi  critical  points.  However,  in  most  electrochemical 
situations  the  electrode  surface  is  polarized  either  positively  or 

negatively,  which  means  that  either  the  cation  or  the  anion  is  strongly 
repelled  from  the  surface,  and  therefore  we  need  to  consider  the 

adsorption  of  either  a  or  b  and  the  solvent  n  on  the  electrode.  This 
implies  a  drastic  simplification  in  the  model,  because  now  we  can 
discuss  at  least  the  case  of  commensurate  adsorption  in  terms  of 
models  that  have  been  solved  analytically,  such  as  the  spin  S=l/2  Ising 
model  and  the  hard  hexagon  gas  model.  The  phase  transitions  predicted 
for  these  models  seem  to  be  reasonable  in  terms  of  the  currently 

available  experimental  evidence.  Consider 


all  sites 
on  the  surf  see 

Here  the  sum  has  two  kind  of  factors,  the  fugacity  of  the 
adsorbed  molecules 
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and  another  factor  which  corresponds  to  the  interaction  of  the 
adsorbed  particles 


g°<r1...rs)=exp[-p[wO(r,...rs) 


(131) 


Equation  (114)  is  the  grand  canonical  partition  function  for  a 
lattice  gas  with  arbitrary  interactions.  In  the  most  general  case,  the 
adsorbate  could  occupy  one  or  more  than  one  single  adsorption  site,  and 
any  adsorbed  molecule  could  interact  with  an  arbitrary  number  of 
neighbors.  This  problem  is  however  quite  untractable,  and  therefore  not 
very  useful.  The  first  approximation  that  comes  to  mind  is  the  Kirkwood 


superposition  approximation 


-rs>  ■  n  s2<'V,j> 


(132) 


Because  of  geometrical  considerations,  this  approximation  is  probably  a 
good  one  for  molecules  with  short  range  interactions  such  as  hard 
spheres,  and  not  a  very  good  one  for  unscreened  charged  panicles.  In 
terms  of  the  effective  potentials,  we  have 


ws(rrrs>-  l  w2<ri'rj> 


1<  i  u<* 


(133) 


With  this  approximation  the  problem  can  be  mapped  onto  the  lattice 
gas  problem  with  arbitrary  interactions,  which  is  still  a  difficult 
problem.  .  If  we  restrict  the  interactions  to  nearest  neighbors  only, 
then  not  only  is  the  problem  a  tractable  one,  but  there  is  a  rather 
extensive  literature  on  cases  that  are  of  physical  interest,  for  which 
the  phase  diagram  of  the  two  dimensional  lattice  gas  is  known,  and 
hence  the  adsorption  isotherm  can  be  deduced. 

Consider  the  case  of  a  simple  salt  dissolved  in  water,  near  a 

metallic  electrode:  in  the  SSM  there  will  be  three  components:  the 
anion,  the  cation  and  the  solvent,  and  the  lattice  atoms.  In  the 

limiting  case  of  the  SSM,  the  sizes  of  the  different  species  play  a 

crucial  role  in  the  possible  ordering  of  the  ad  layers  at  the 


interface.  It  will  be  convenient  to  picture  the  tons  as  having  a  hard 
sphere  core  with  a  diameter  oa>  and  the  solvent  as  having  a  hard 
core  with  diameter  cn.  The  lattice  spacing  of  the  metal  surface  is  d, 

and  because  it  is  the  most  stable  surface,  we  will  restrict  ourselves 
to  the  (111)  surface  of  the  fee  crystals,  or  the  (100)  face  of  the  hep 
crystal,  that  is,  the  triangular  lattice.  In  the  most  general  case,  all 
three  components  can  be  adsorbed  competitively,  and  this  situation  can 
give  rise  to  very  complex  phase  diagrams.  In  most  electrochemical 

situations  the  electrode  surface  is  polarized  either  positively  or 
negatively,  which  means  that  either  the  cation  or  the  anion  is  strongly 
repelled  from  the  surface,  and  therefore  we  need  to  consider  the 
adsorption  of  either  a  or  b  and  the  solvent  n  on  the  electrode.  This 

implies  a  drastic  simplification  in  the  model,  because  now  we  can 
discuss  at  least  the  case  of  commensurate  adsorption  in  terms  of 
models  that  have  been  solved  analytically,  such  as  the  spin  s=l/2  Ising 
model  and  the  hard  hexagon  gas  model.  The  phase  transitions  predicted 

for  these  models  seem  to  be  reasonable  in  terms  of  the  currently 
available  experimental  evidence. 

There  are  two  cases  of  chemisorption  of  electrochemical  interest: 
In  the  first  one,  the  charge  of  the  adsorbate  is  neutralized  by  the 
electrons  in  the  metal,  and  this  means  that  the  interactions  between 
neighbors  on  the  surface  is  attractive.  If  we  are  far  from  the  point  of 
zero  charge,  and  the  metal  is  negatively  charged,  the  contact 
probability  of  the  anions  a  is  zero  for  all  practical  purposes  and  we 
have  only  the  cation  b  or  the  solvent  n  on  the  adsorption  sites.  The 
problem  is  then  reduced  to  a  spin  s=l/2  Ising  model  with  ferromagnetic 
interactions.  In  the  second  case  there  is  no  discharge  of  the  adions  by 
the  metal.  The  interactions  between  the  ions  of  the  same  sign  is 
clearly  repulsive,  so  that  the  nearest  neighbor  sites  to  an  occupied 
site  are  not  going  to  be  occupied. 

The  possible  existence  of  phase  transitions  was  discussed  in  the 
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work  of  Huckaby  and  Blum  .  The  observation  is  that  the  SSM  model  maps 
the  three  dimensional  interface  onto  a  two  dimensional  lattice  problem. 
The  phase  behavior  in  the  interphase  is  determined  by  the  mapping  of 
the  parameters,  and  exact  conditions  on  the  existence  of  phase 
transitions  can  be  given. 

The  simplest  electrolyte  has  three  components  that  can  be  adsorbed 


35 


onto  the  sites  of  the  metal:  the  anion,  the  cation  and  the  solvent.  The 
possible  phase  behavior  of  such  a  system  has  been  discussed  in  the 
literature  and  although  there  are  some  exact  results  on  the  phase 
behavior  of  the  adsorbed  layer,  the  study  of  this  system  remains  an 
open  problem.  The  simpler  case  of  two  component  adsorption  is  discussed 
in  more  detail. 


THE  THREE  STATE  ADSORPTION  MODEL 

This  is  the  case  of  a  simple  salt  dissolved  in  a  solvent  like 
water.  We  call  the  ions  a  and  b,  and  the  solvent  n.  since  this  is  a 
dense  system,  then  the  sdcky  sites  will  never  be  empty,  they  will  be 

occupied  by  either  by  a,  b  or  n.  Now  there  are  various  possibilities, 

because  the  size  of  adsorbate  may  be  bigger  or  smaller  that  the 
site- site  separation.  In  the  first  case  the  adsorbate  will  exclude  not 

only  us  own  site  but  also  the  neighboring  ones,  and  the  problem 

becomes  the  hard  hexagon  problem  of  Baxter  when  the  interactions  of  the 
neighboring  sites  are  ignored  altogether.  If  the  adsorbate  is  smaller 
than  the  lattice  site  separation  then  it  can  been  shown  that  the 
problem  can  be  reduced  to  the  two  component  lattice  gas:  This  case  was 
recently  studied  by  Rikvold  .  In  fact  we  have,  by  direct  comparison  of 
the  partition  function  and  the  hamiltonian  of  the  three  state  lattice 
gas 


(134) 

where  we  have  used  Rikvold’s  notation:  the  operators  c-  are  1  when  sue 
i  is  occupied  by  x,  and  0  otherwise.  In  our  case  x  is  either  a  or  b, 
and  the  solvent  n  counts  as  the  empty  site.  The  interactions  are 
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(135) 


♦a.=  waa  ’  wnn  =  ln  (  4.<Vrj>/gnn(r/j)  ) 

i  ,j  are  neighboring 
sites 

and  similar  relations  for  0^  and  0^-  The  last  two  sums  in  (134) 
correspond  to  the  single  occupation  of  a  site  by  x  : 

Pa=(l/(3)  ln  *a/*n=  (1/(3)  ln 

(136) 

and  a  similar  expression  for  (i^. 

The  possible  arrangements  on  the  triangular  lattice  of  a  three 

component  mixture  is  quite  complicated  and  has  been  studied  recently  by 
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Collins  et  al.  (see  also  Y.Saito  ).  There  are  10  different  ordered 
phases  on  the  surface,  when  only  next  nearest  neighbors  are 
interacting.  Of  course  the  picture  is  even  more  complicated  for  longer 
ranged  interactions.  The  geometrical  arrangement  is  either  (lxl),  and 
there  are  three  phases  of  this  kind  which  correspond  to  pure  a,  pure  b 
and  pure  solvent  n.  These  are  dense  phases  and  correspond  to  total 
coverage  (0=1)  for  each  of  these  components.  These  are  limiting  cases 
when  there  is  very  strong  adsorption  of  any  of  the  components  and 
attractive  interaction  between  atoms  of  the  same  kind.  If  we  consider 
ions  with  their  charge,  then  because  of  the  Coulombic  repulsion  the 
interaction  is  repulsive  or  antiferromagnetic  and  the  dense  phases  with 
either  pure  a  or  pure  b  are  to  be  excluded.  In  the  case  of  the 
underpotential  deposited  dense  monolayers.  This  is  explained  by  the 
fact  that  the  metal  cations  are  not  charged  at  the  surface,  and 
therefore  do  not  repell  each  other. 

The  remaining  7  phases  are  mapped  on  the  (43x43)  lattice.  The  unit 
cell  of  this  lattice  contains  three  sites  of  the  original  triangular 
lattice,  and  appears  rotated  by  30*  with  respect  to  it.  Each  of  the 
three  sites  can  be  occupied  by  either  of  the  components  of  the  mixture, 
but  because  of  the  symmetry  (geometrical  degeneracy)  of  the  phases,  the 
number  of  distinct  coverings  is  reduced  from  24  to  7.  In  table  1,  which 


r  -4-'- 1 
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is  taken  directly  from  Collins’s  work  we  show  the  coverage,  energy  and 
degeneracy  of  each  of  the  ordered  phases. 

There  are  two  cases  chemisorption  of  electrochemical  interest:  In 
the  first  one,  the  charge  of  the  adsorbate  is  neutralized  by  the 
electrons  in  the  metal,  and  this  means  that  the  interactions  between 
neighbors  on  the  surface  is  attractive.  If  we  are  far  from  the  point  of 
zero  charge,  and  the  metal  is  negatively  charged,  the  contact 
probability  of  the  anions  a  is  zero  for  all  practical  purposes  and  we 
have  only  the  cation  b  or  the  solvent  n  on  the  adsorption  sites.  The 
problem  is  then  reduced  to  an  spin  s=  1/2  Ising  model  with  ferromagnetic 
interactions,  which  will  be  discussed  in  the  next  section. 

In  the  second  case  there  is  no  discharge  of  the  adions  by  the 
metal.  The  interactions  between  the  ions  of  the  same  sign  is  clearly 
repulsive,  and  that  of  opposite  sign  is  attractive:  we  have  then 

fab<0'an<i  W50  (137) 

This  is  definitely  the  antiferromagnetic  case.  Let  us  assume 
for  simplicity  that  we  are  dealing  with  the  symmetric  case  in  which 


♦aa  =  *1 


(138) 


Then  there  are  two  possible  situations:  either  the 
repulsive  interactions  predominate  and 


'  ^aa  >  ^ab  (casc  l) 


or  the  ab  attractions  predominate 


(139) 


-  <>aa  4  ♦ab  (case  n> 

(140) 

The  phases  in  case  I  for  a  possible  set  of  parameters  is  pictured 
in  figure  7.  Here  we  may  have  phases  ranging  from  pure  a,  for  very 
large  negative  pa  to  pure  b  for  very  large  negative  (i^.  In  the  example 
of  case  II  the  ab  attraction  predominates  so  as  to  render  the  dilute 
phases  with  9j?=l//3,  2/3  thermodynamically  unstable.  This  case  is 
pictured  in  figure  8. 

Complicated  phase  diagrams  may  exist  in  the  neighborhood  of 
the  point  of  zero  charge  for  systems  in  which  there  is  strong 
chemisorption  of  some  of  the  species.  In  the  absence  of  specific 
adsorption  one  would  expect  that  either  p or  p^  is  large  ,  but  not 
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both  at  the  same  time.  In  these  cases  we  are  in  the  asymptotic  regime 
where  the  system  can  be  created  by  two  state  models.  This  will  be  done 
in  the  next  section. 

THE  TWO  STATE  ADSORPTION  MODEL 
We  specialize  now  to  the  case  in  which  only  one  of  the  ionic 
species  of  the  solution  or  the  solvent  is  adsorbed.  In  this  case  the 
structure  of  the  phase  diagrams  is  simpler,  which  at  the  same  time 
allows  for  a  more  general  discussion,  because  we  can  include  the  cases 
in  which  the  size  of  the  adsorbate  is  variable. 

The  following  discussion  is  taken  from  the  work  of  Huckaby  and 
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Blum  .  We  consider  a  system  in  which  we  have  a  triangular  lattice  of 
sticky  sites  on  a  hard  plane  surface.  The  spacing  between  the  sites  of 
the  lattice  is  d.  This  surface  is  in  contact  with  a  solution.  Only  two 
states  of  occupation  are  allowed:  the  sites  are  either  occupied  by  an 

ion  or  by  the  solvent,  or  alternatively  in  the  case  of  a  pure  fluid,  by 
a  fluid  particle  or  none.  The  fluid  particles  have  an  exclusion 
diameter  a,  which  may  or  may  not  be  associated  with  a  hard  core 

potential.  Otherwise  the  interactions  are  arbitrary.  We  assume  however 
that  the  pair  correlations  on  the  surface  decay  sufficiently  fast  so 
that  we  need  to  take  into  account  first  neighbor  interactions  only. 

There  are  two  possible  situations:  If  the  adsorbate  diameter  o  is 
smaller  than  the  lattice  spacing  d,  then  there  are  two  possible  phases, 
a  dense,  crystalline  one  and  a  dilute  disordered  one.  Three  is  a  first 

order  transition  between  them.  If  the  adsorbed  particles  exclude  all 
next  nearest  neighbors,  that  is  when 

d<c<43d  (141) 
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then  the  problem  is  exactly  analogous  to  the  hard  hexagon  problem  of 
Baxter.  In  this  case  there  is  a  second  order  phase  transition  between 
an  ordered  43x43  phase  and  a  disordered  one. 

SINGLE  SITE  OCCUPANCY 

Consider  again  equation  (114)  together  with  the  superposition 
approximation  (133): 
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(142) 


*sn  p  <ri> 
Z/Z°=  Y  - Hi - 

sj 

All  site* 
on  the  surface 


where  the  subscript  nn  means  that  the  pairs  (r^,ij)  *n  the  sum  inside 
the  square  brackets  are  nearest  neighbors.  As  was  done  before  (see 
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eq.134),  we  can  write  (142)  as  a  lattice  gas  partition  function 


Z/Z°  -  Y.  exp['P*Z‘i  'j  +  ‘i  ] 


^=0,1 


(143) 


where  t-  is  the  occupation  number  of  site  i  which  can  be  either  0  for 
the  site  when  it  is  occupied  by  the  solvent  (or  empty)  and  1  when  it  is 
occupied.  Furthermore 


Pw  =  -  In  g^(d) 
PH  =  ln[Xp°(0)] 


(144) 


(145) 


This  partition  function  can  be  mapped  onto  an  Ising  model  with 
spin  variables  s.=±l  by  means  of  the  transformation 

si=2ti- 1  or  tpfSj+D/2  (146 


(146) 


In  this  case  (143)  becomes  the  partition  function  for  the  Ising 


model 


Z/Z°  =  C  £  exp [-§21 1  (siSj  -2sf2Sj)  +4  ] 


si=±  1 


(147) 


where  C  is  a  constant  that  is  irrelevant  to  our  calculation.  w/4=J  is 
the  Ising  parameter,  which  is  bigger  than  zero  in  the  ferromagnetic 
case,  and  smaller  than  zero  in  the  antiferromagnetic  case.  The  variable 
that  plays  the  role  of  the  external  magnetic  field  is  h,  given  by 
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2h  =  p-wq/2  (148) 

where  q  is  the  coordination  number  of  the  lattice.  In  the  case  of  the 
triangular  lattice  q=6.  The  two  state  ferromagnetic  Ising  model  has  a 
first  order  phase  transition  when  h=0  .  From  (148)  this  means 


p=wq/2  (149) 

and  using  the  definitions  (144)  and  (145)  we  get  the  exact  condition  for 
phase  transitions  when  w>0 

?=Xp°(0)=[g°(d)r3  (150) 

The  preceding  analysis  can  be  illustrated  by  a  fluid  of  hard 
spheres  in  contact  with  a  sticky  triangular  lattice  of  spacings  d=G.  In 
this  case  a  good  estimate  of  both  the  contact  density  and  the  pair 
distribution  function  are  obtained  from  the  Percus  Yevick  theory 


p°(0)=p  mini 

(1-Tir 


(151) 


(l+Ti/2) 

(i-n)2 


(152) 


Ti=(l/6)Jtpo3  (153) 

where  T|  is  the  fraction  of  occupied  volume.  Replacing  into  (150)  we  get 

3  o-n)8 

X  =  (l/6)lt<r - r  (154) 

T]  (1+2T1)  (l+n/2) 

This  relation  is  a  necessary  but  not  sufficient  condition  for  the 
occurrence  of  phase  transitions. 

A  sufficient  condition  for  the  occurrence  ot  a  phase  transition 

can  be  obtained  from  the  classic  work  of  Potts  on  the  magnetization  of 

95 

the  ferromagnetic  Ising  model  on  the  triangular  lattice  .  The 
magnetization  in  the  Ising  model  is  the  difference  in  the  number  of 
spins  up  pu  and  spins  down  p^ 

I(x)=pu-Pd  (155) 


or  using  (146)  I(x)=  20-1  (156) 

where  0  is  the  fraction  of  occupied  sites.  Following  Potts  we  use  the 
variable 


x.  e-Pw/2= 


(157) 
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The  result  for  the  spontaneous  magnetization  is 


I(x)=  [l- 


(l+3xz)(l-x): 


(158) 


Using  (156)  and  (157)  yields  the  fraction  of  occupied  sites  as  a 
function  of  the  pair  interaction  function 


0  =  1/2  -(1/2)  1 


16  g^(d) 

(g^(d)+3)(  g^(d)-l)- 


(159) 


The  critical  value  of  0  occurs  when  0=1/2.  Solving  (159)  for  the 
contact  pair  correlation  function  yields 


g°(d)=3 

From  (150)  we  get  the  value  for  the  critical  sticky  parameter  A. 


(160) 


A.  p°(0)=l/27  (161) 

In  our  example  of  hard  core  fluid  in  the  Percus  Yevick 
approximation,  we  get  using  these  equations  together  with  (151)  and 
(152)  we  get  the  critical  value  of  the  excluded  volume  fraction 

Tlcrir0-3712  062) 

and  the  sticky  parameter 

X^.01185*3  (163) 

This  system  undergoes  a  first  order  phase  transition.  The 
isotherms  for  various  values  of  the  parameter  A,  are  shown  in  figure  9. 
These  exact  results  are  in  qualitative  agreement  with  the  mean  field 

O  A 

theory  of  Badiali  et  al:  ,  where  the  first  order  phase  transition  is 
also  predicted.  But  the  quantitative  agreement  is  not  good,  which 
illustrates  the  pitfalls  of  mean  field  theory.  In  the  computer 
simulations  of  Caillol  et.  al.^  the  conditions  for  the  occurrence  of  a 
first  order  phase  transition  are  not  met,  because  the  adsorption  sites 
are  of  finite  size,  and  for  that  reason  the  occupancy  of  a  site  may 
prevent  nearest  neighbor  occupation.  In  this  case  we  expect  a  second 
order  phase  transition  to  occur. 

When  there  are  longer  ranged  interactions,  beyond  the  nearest 

07 

neighbor  interactions,  Dobrushin  has  shown  that  the  first  order  phase 
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transition  still  occurs.  However,  an  exact  relation  such  as  (150)  is  not 
available. 

As  was  mentioned  already  in  the  electrochemical  case  the  contact 

pair  correlation  function  of  ions  of  equal  sign  is  practically  zero, 

because  of  the  Coulomb  repulsion  which  prevents  ions  of  equal  sign  to 

approach  each  other.  However  condensed  phases  in  the  adlayers  are 

observed  in  electrochemistry.  In  particular  the  underpotential 

deposition  of  some  metals  on  electrodes  occurs  at  certain  very  well 

98 

defined  values  of  the  potential  bias  .  For  example  the  deposition  of  Cu 

on  the  Au  (111)  face  forms  two  phases  according  to  the  deposition 

99  100  101 

potential  .  These  phases  have  been  observed  ex-situ  and  in  situ  .  At  a 

lower  potential  a  dilute  ordered  <13x4  3  phase  is  formed.  At  a  higher 

potential  a  dense  commensurate  phase  is  formed.(  see  figure  10).  It  is 

clear  from  the  above  considerations  that  in  the  dense  adlayer  case  the 

ions  must  be  discharged,  because  then  they  would  form  a  metallic  bond, 

which  makes  w  positive  ,  and  therefore  ferromagnetic.  This  is 

supported  by  the  features  of  the  EXAFS  spectra.  In  the  high  density 

phase  the  near  edge  structure  corresponds  to  that  of  metallic  copper, 

which  has  a  characteristic  double  peak  (figure  11).  The  dilute  4  3x4  3 

phase  has  the  white  line  characteristic  of  the  charged  ions.  We  may 

assume  then  that  in  this  case  the  Cu  retains  pan  of  its  charge,  so 

that  the  interactions  are  in  this  case  repulsive,  which  corresponds  to 

the  antiferromagnetic  case. 

MULTIPLE  SITE  OCCUPANCY:  THE  HARD  HEXAGON  CASE 

If  the  occupation  of  one  site  in  the  triangular  lattice  also 

excludes  the  nearest  neighbors,  then  the  problem  is  equivalent  to  the 
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hard  hexagon  problem  of  Baxter  .  This  problem  can  be  solved  when  the 

interactions  between  the  hard  hexagons  are  neglected. 

The  thermodynamics  of  the  hard  hexagon  model  was  recently  worked 
102 

out  by  Joyce  .  In  terms  of  the  lattice  fugacity 

*=Xp°(0)  (130) 

The  system  undergoes  a  second  order  phase  transition  between  an 
ordered  solidlike  phase  and  a  disordered  one.  The  transition  occurs 
when 

*  =(1 1+54  5)/2= 11.09  (164) 

Wall 
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The  fraction  of  occupied  sites  is 


0crit  =(5-J5)/10=0.2764  (165) 

We  remark  that  in  the  limit  of  highest  possible  density  the 
occupied  sites  fraction  is  0=1/3. 


For  the  low  density  phase  the  equation  of  state  is 
1  r  2  1/2  /  1/2 

+  Q2  -  %[  2Q3  +  2<32  <*1 


Qq=  1-58+50 2 
Qp  Qq(  1-0+02) 


Q2=(  1  -  20)  ( 1  - 1 10+4402  -7703+6604-  3305+l  196) 


Q3=l  -  160+lO602-37803+8O304-lO8O05+96206 
-  5  7  607 +2 1 90  8  -  5O09+ 1 00 1 0 


while  for  the  ordered  phase 

2(2 - 30) ( 1-0)3 

?=  - - 

|-l  +  120-4502+6603-3304-  [(l-50+502)3(l-90+902)j  j 

Consider  the  case  of  an  electrolyte  in  the  neighborhood  of  a 
charged  electrode.  The  simplest  model  is  the  restricted  primitive  model 
of  a  continuum  dielectric  and  a  smooth  imageless  surface.  The  contact 
probability  is  computed  from  the  Gouy-Chapman  theory  discussed  in 
section  I. 

g2(d)=p.(0)/p=  exp(-z.y.(0)]  (168) 

where  we  are  using  (7)  together  with  the  definition 

\)T(0)=  pep.(0)  (169) 

we  get 

?=Xp.exp[-z.\j/.(0)I  ( 1 70) 

where  yXO)  is  the  reduced  potential  at  the  inner  Helmholtz  plane. 
Equating  (170)  to  either  (166)  or  (167)  produces  a  suggestive  form  of 
an  adsorption  isotherm  that  is  shown  in  figure  12. 
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(166) 


(167) 


From  (164)  and  (170)  we  deduce  that  in  reduced  units  there  is  a 
critical  potential  at  which  a  phase  transition  occurs  between  an 
ordered  and  a  disordered  phase 

z  y  (0)=ln[Xp  j-2.406  (l71 

i  cm  i 

At  this  point  a  second  order  phase  transition  occurs  in  the 
surface,  which  is  also  seen  in  the  adsorption  isotherm. 
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TABLE  1 


Configuration 

6a 

9b 

0n 

Energy 

Degeneracy 

Dense  Phases  (lxl) 

a 

a  a 

1 

0 

0 

3<t>aa+^a 

1 

b 

b  b 

0 

1 

0 

^bb+^b 

1 

n 

n  n 

0 

0 

1 

0 

1 

Diluted  Phases  (>I3x>I3) 

b 

a  a 

2/3 

1/3 

0 

*ab+(i)*bb+  (i)<*>aa 

♦<i>  V(!> 

3 

a 

b  b 

1/3 

2/3 

0 

^ab+{  i  )<^aa+  (i}<l>bb 
*(\) 

3 

n 

a  a 

2/3 

0 

1/3 

^aa  +  ^  3  ^  ^a 

3 

n 

b  b 

0 

2/3 

1/3 

*bb  ♦<!>  ^b 

3 

a 

n  n 

1/3 

0 

2/3 

<!>  ^a 

3 

b 

n  n 

0 

1/3 

2/3 

3 

a 

b  n 

1/3 

1/3 

1/3 

^ab  +  V(S}  ^b 

6 
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FIGURE  CAPTIONS 


Figure  1:  The  primitive  model  of  the  electric  double  layer.  The  metal 
side  is  represented  by  a  smooth  hard  wall.  The  ions  are  charged  hard 
spheres.  The  solvent  is  a  dielectric  continuum. 

★ 

Figure  2:  Density  profile  for  the  test  case  with  o  =0.7 

*  , 

(ref.34,c  =a‘E/e).  The  line  is  the  result  of  Ballone,  Pastore  and  Tosi, 
ref.  54. 

Figure  3:  same  as  fig.2.  Results  from  Nielaba  and  Forstmann  ,  ref.56. 

Figure  4:  same  as  fig.2.  Results  from  Caccamo,  Pizzimenti  and  Blum, 

ref.59. 

Figure  5:  same  as  fig.2.  Results  from  Plischke  and  Henderson,  ref.66. 

Figure  6:  same  as  fig.2.  Results  from  Kjellander  and  Marcelja,  ref.  68. 

Figure  7:  Phase  diagram  for  the  three  state  adsorption  model  case  I 
(from  Rikvold  et  al.  ref.  89). 

Figure  8:  same  as  figure  7.  case  II. 

Figure  9:  Adsrption  isotherms  for  hard  spheres  in  the  mean  field 
approximation  for  a  triangular  lattice  of  sticky  sites,  (ref.84). 

Figure  10:  Cyclic  potential  curves  for  underpotential  deposition  of  Cu 
on  Au(lll).  The  first  peak  (low  density  phase),  corresponds  to  a 
coverage  0=0.3.  The  second  peak  corresponds  to  a  dense,  commensurate 
adlayer  (ref.99). 

Figure  11:  Comparison  of  the  near  edge  adsorption  peaks  for  the 
underpotential  deposited  Cu  on  Au  (111),  (a)  spectra  of  the  low  density 
phase.  The  near  edge  structure  has  a  single  high  peak,  similar  to  that 
of  CuSO^,  shown  in  (b).  The  high  density  phase  spectrum  (c)  has  a 
double  peak,  similar  to  that  of  Cu  foil  (d). 

Figure  12:  The  fraction  of  occupied  sites  0  as  a  function  of  the 
applied  potential.  The  units  are  adimensional.  There  is  a  second  order 
phase  transition  at  0=0.2764. 
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